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Abstract 

A search for charged excited leptons decaying into a lepton and a photon has been performed 
using approximately 680 pb _1 of e + e~ collision data collected by the OPAL detector at LEP 
at centre-of-mass energies between 183 GeV and 209 GeV. No evidence for their existence was 
found. Upper limits on the product of the cross-section and the branching fraction are inferred. 
Using results from the search for singly produced excited leptons, upper limits on the ratio 
of the excited lepton coupling constant to the compositeness scale are calculated. From pair 
production searches, 95% confidence level lower limits on the masses of excited electrons, muons 
and taus are determined to be 103.2 GeV. 
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1 Introduction 



Models in which fermions have substructure attempt to explain, among other things, the well- 
ordered pattern of fermion generations observed in nature. The existence of excited states of the 
Standard Model fermions would be a natural consequence of fermion compositeness. Excited 
leptons could be produced in e + e~ collisions and are expected to decay via the emission of a 
gauge boson (7, Z° or W^) ||]. 

This paper presents results from a search for excited electrons (e*), muons (//*) and tau 
leptons (t*) decaying electromagnetically, using data collected by the OPAL experiment at 
LEP. Results presented in this paper are obtained using a larger data sample, as well as a 
significantly improved analysis, compared to our previous publications ||. Searches for excited 
leptons have also been performed by other LEP collaborations || and by the HERA experiments 
in electron-proton collisions @. 

At LEP, excited leptons could be produced in pairs or in association with a Standard Model 
lepton. Both single and pair production of excited leptons proceed through s-channel photon and 
Z° diagrams. In addition, ^-channel photon and Z° exchange diagrams contribute to the single 
and pair production of excited electrons. The ^-channel contribution, although expected to be 
negligible for pair production, causes excited electrons to be singly produced predominantly in 
the forward region with the recoiling electron outside the detector acceptance. Thus, in addition 
to final states containing two leptons and one or two photons (££j, ££jj), a separate search for 
the single production of excited electrons with an undetected electron (e7) is also performed. 

The results presented in this paper are interpreted in the context of the framework described 
in |],||]. In this phenomenological model, the interaction between excited leptons and a gauge 
boson (£*£*V), which largely determines the cross-section for pair production of excited leptons, 
is vector-like. The single production cross-section and branching fractions of excited leptons 
are determined by the strength of the ££*Y coupling. This interaction can be described by the 
following SU(2)xU(l) gauge invariant effective Lagrangian 0,0 



2A 



£l + hermitian conjugate, 



where a^ u is the covariant bilinear tensor, r denotes the Pauli matrices, Y is the weak hyper- 
charge, and B^ v represent the Standard Model gauge field tensors and the couplings g, g' 
are the SU(2) and U(l) coupling constants of the Standard Model. The compositeness scale is 
set by the parameter A which has units of energy. Finally, the strength of the ££*\ coupling 
is governed by the constants / and /'. These constants can be interpreted as weight factors 
associated with the different gauge groups. The values of / and /' dictate the relative branching 
fractions of excited leptons to each gauge boson. The branching fraction of electromagnetically 
decaying excited charged leptons is significant for most values of / and /' except in the case 
where f = —f which entirely forbids this particular decay. As a result of the clean characteris- 
tic signatures expected, the photon decay constitutes one of the most sensitive channels for the 
search for excited leptons, even for values of / and /' for which other decay modes dominate. 
To reduce the number of free parameters it is customary to assume either a relation between 
/ and /' or set one coupling to zero. For easy comparison with previously published results, 
limits calculated in this paper correspond to the coupling choice / = /'. This assignment is a 
natural choice which forbids excited neutrinos from decaying electromagnetically. For this par- 
ticular coupling choice, the electromagnetic branching fraction of charged excited leptons drops 
smoothly from 100% for masses below the Z° and W ± mass thresholds to about 30% for masses 
in excess of 200 GeV. 



3 



bin range 
(GeVJ 


< V s > 
(GeVJ 


L 

(pb ) 


178.00 - 186.00 


182.7 


63.8 


186.00 - 190.40 


188.6 


183.2 


190.40 - 194.00 


191.6 


29.3 


194.00 - 198.00 


195.5 


76.5 


198.00 - 201.00 


199.5 


76.9 


201.00 - 203.75 


201.9 


44.5 


203.75 - 204.25 


203.9 


1.5 


204.25 - 204.75 


204.6 


9.7 


204.75 - 205.25 


205.1 


60.0 


205.25 - 205.75 


205.4 


3.6 


205.75 - 206.25 


206.1 


14.3 


206.25 - 206.75 


206.5 


1U 1 .6 


206.75 - 207.25 


206.9 


5.7 


207.25 - 207.75 


207.5 


0.5 


207.75 - 208.25 


208.0 


7.2 


> 208.25 


208.3 


0.5 




684.4 



Table 1: The luminosity weighted mean centre-of-mass energy and integrated luminosity of 
each energy bin. 



2 Data and Simulated Event Samples 

The data analysed were collected by the OPAL detector at centre-of-mass energies ranging 
from 183 GeV to 209 GeV during the LEP runs in the years 1997 to 2000. The search for excited 
leptons is based on a total of 684.4 pb -1 of data for which all relevant detector components were 
fully operational. For the purpose of accurately interpreting the results in terms of limits on 
excited lepton masses and couplings, the data are divided into 16 centre-of-mass energy bins 
analysed separately. The energy range, luminosity weighted mean centre-of-mass energy and 
integrated luminosity of each bin are summarised in Table |l|. The uncertainty on the measured 
beam energy is approximately 25 MeV § and is correlated between centre-of-mass energy bins. 
In addition to the high energy data, approximately 10 pb -1 of calibration data collected in 1997- 
2000 at a centre-of-mass energy near the Z° mass were used to study the detector response. 

Distributions of kinematic variables and selection efficiencies for excited leptons were mod- 
elled using samples of simulated events obtained using the EXOTIC || Monte Carlo event 
generator. The matrix elements |||l0| implemented in EXOTIC include all the spin correlations 
in the production and decay of excited leptons. 

The Standard Model processes at different centre-of-mass energies were simulated using 
a variety of Monte Carlo event generators. Bhabha events were simulated using the BH- 



WIDE JllJ and TEEGG [|12[] generators, muon and tau pair events using both KORALZ [|l 



and KK2F |]l4f, e + e — ► qq(7) events using PYTHIA |Tq| and KK2F, four fermion processes 



using KORALW [16| and grc4f fi7|| , di-photon production using RADCOR [|lq| , and two-photon 



events (e+e" e+e"77 e+e"//) using VERMASEREN gg|, PHOJET @ and HER- 
WIG [^TJ. Each of the simulated event samples was processed through the OPAL detector 
simulation program p2| and analysed in the same way as data. 
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3 Event Selection and Kinematic Fits 



Events are reconstructed from tracks and energy clusters defined by requirements similar to 



those described in [23|. The background from multihadronic events is substantially reduced by 



requiring at least one but no more than six tracks in an event. Furthermore, the ratio of the 



number of good tracks, as defined in [23|, to the total number of tracks reconstructed in the 
detector must be greater than 0.2 in order to reduce background from beam-gas and beam-wall 
collisions. Cosmic ray events are suppressed using information from the time-of-flight counters 
and the central tracking chamber pj]. 



Tracks and energy clusters in an event are grouped into jets using a cone algorithm [25| with 
a cone half-angle of 0.25 radians and minimum jet energy of 2.5 GeV. The parameters defining 
a jet were chosen to maximise the signal efficiency over the broadest possible range of excited 
lepton masses and centre-of-mass energies. Events are required to contain between two and four 
jets. Jets are classified as leptons or photons using the criteria described below applied in the 
same order as given in the text. 

Photon candidates must have a minimum energy deposited in the electromagnetic calorimeter 
equivalent to 5% of the beam energy. A photon jet must either contain no tracks or be identified 



as a photon conversion using a neural network technique [26]. Jets in which the most energetic 
track has a neural network output greater than 0.9 and the energy deposited in the hadronic 
calorimeter is less than 10% of the beam energy are defined to be photon conversions. All photon 
candidates must lie within | cos 9\ < 0.9 to avoid poorly modelled regions of the detector []. The 
energy and direction of each photon candidate are determined from the energy and position of 
the energy cluster in the electromagnetic calorimeter. 

Muon candidates are jets containing exactly one track with associated hits in the muon 
detectors or hits in the hadronic calorimeter consistent with the particle being a muon [ p^| . 
Muons, unlike electrons and photons, do not shower while traversing additional material present 
in the forward region of the detector. Muon candidates are thus allowed to lie within a larger 
angular acceptance of | cos 6\ < 0.95. The direction of each muon candidate is given by the polar 
and azimuthal angles of the track and the momentum is calculated from the track curvature and 
polar angle. 

A jet is identified as an electron if it contains exactly one track satisfying one of the following 
two requirements: the ratio of the electromagnetic energy to the track momentum (E/p) lies 
between 0.8 and 1.4 or the track has an output greater than 0.9 from a neural network developed 
to identify electrons [26]. Electron candidates are also required to lie within |cos#| < 0.9 to 



avoid poorly modelled regions of the detector. The energy of each electron candidate is taken 
to be the energy deposited in the electromagnetic calorimeter while the direction is given by the 
polar and azimuthal angles of the track. 

Finally, unidentified jets containing at least one track and lying within | cos 6*| < 0.95 are 
considered to be tau candidates. Jets in the region 0.90 < | cos#| < 0.95 which would satisfy the 
electron or photon requirements are discarded from the sample of tau candidates. Jets identified 
as tau candidates are mostly hadronically decaying taus. Tau leptons decaying leptonically are 
tagged as electrons or muons by the criteria described above. The polar and azimuthal angles 
of tau candidates are given by the axis of the jet, corrected for double-counting of tracks and 



energy clusters |28j. 



The different selections used to identify the final states of interest are described in the 
following sections and the results are summarised in Table ||. 



lr The OPAL coordinate system is denned to be right-handed, with the z-axis pointing along the electron beam 
direction and the x-axis pointing toward the centre of the LEP ring. Thus the polar angle 9 used in this paper 
refers to the angle with respect to the electron beam direction and the azimuthal angle cf>, the angle measured 
with respect to the x-axis. 
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3.1 Selection of ££77 Final States 



Events containing two lepton candidates of the same flavour and two identified photons are 
considered as candidate events for the pair production of excited leptons. In addition, events 
containing two leptons of different flavours and two photons are considered as excited tau can- 
didates. In order to reduce the background from Standard Model processes, additional selection 
criteria are applied to the different types of candidate events. 

The quantity R v [ s is defined to be the sum of the energy of the particles considered for 
a given event final state, divided by the centre-of-mass energy. This quantity is required to 
exceed 0.8 for ee77 and ^^77 candidates and 0.4 for TT77 candidates. This criterion reduces 
the background from two-photon events. It also decreases the contamination from qq events in 
the TT77 sample. Figure ||(a-c) shows the i? v i S distributions obtained using the entire data set 
for each type of candidate event. The observed discrepancy at small values of i? v i s in the TT77 
sample corresponds to a region where the background is dominated by two-photon events and 
does not affect the analysis as the events of interest lie in a region of i? v ; s that is well modelled. 
A similar mis- modelling is present in distributions of ££j and e7 candidate events. 

The remaining background in the ee77 and /^77 samples comes almost entirely from e + e~ 
and [i + {J-~ events with additional photons. The background in the TT77 sample consists mostly 
of t + t~ events with more than one radiated photon, and a small fraction of qq events. 

3.2 Selection of ££j Final States 

Events containing two lepton candidates of the same flavour and at least one identified photon 
are considered as candidate events for the single production of excited leptons. In addition, 
events with two leptons of different flavours and at least one photon are considered as excited 
tau candidates. If more than one photon is identified in the event, the most energetic photon is 
chosen and the other photon is ignored. Events selected as candidates for the pair production 
of excited leptons are also considered as single production candidates. 

To reduce the background from two-photon events the quantity -R v ; s must be greater than 
0.8 for ee7 and ^7 candidates, and greater than 0.4 for ttj candidates. The R v i s distributions 
of each type of candidate event are shown in Figure |l](d-f ) for data from all the centre-of-mass 
energies combined. 

The dominant e + e~ background in the ee7 final state is reduced by requiring that the angle 
between the most energetic electron and photon (9 e7 ) be greater than 90°. The cos 6 ei distribu- 
tion obtained using data from all centre-of-mass energies combined is shown in Figure ||(a). 

Background from both e + e~ and fi + n~ events in the ttj sample is reduced by requiring the 
total energy deposited in the electromagnetic calorimeter to be between 20% and 80% of the 
centre-of-mass energy. Finally, the polar angle of the missing momentum vector for the particles 
considered in the TT7 final state must lie within | cos# m i ss | < 0.9. This requirement reduces the 
contamination from qq events. Figures |2](b,c) show both the total electromagnetic energy and 
I cos# m i SS | distributions of ttj events before applying each cut. 

After this selection, the remaining background in the ee7 and hhj samples consists almost 
entirely of e + e~ and events with an additional photon. The background in the TT7 sample 

consists mostly of t + t~ events with one radiated photon and a small fraction of radiative e + e~, 
and qq events. 

3.3 Selection of e7 Final State 

A separate selection for events with one electron and one photon was developed to increase the 
efficiency of the search for singly produced excited electrons where one electron travels in the 
forward region outside the detector acceptance. 
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Candidate events are required to contain at least one photon and at least one electron 
candidate. Additional jets, if present, are ignored. Since the ee7 and e7 final states are combined 
to calculate a limit on the single production of excited electrons, it is important to ensure that 
events are not double-counted. All events that are selected by the set of general requirements 
discussed at the beginning of Section ||, but that fail the ee7 selection are considered as possible 
e7 candidates. 

To reduce the two-photon background, the quantity i? v i s must satisfy i2 v ; s > 0.4. The angle 
between the electron and photon {Q e -y) is also required to be greater than 90°. Further reduction 
of the dominant e + e~ background is achieved by requiring that the measured polar angle of 
the photon satisfies | cos# 7 | < 0.8 and by rejecting events where the photon is identified as a 
conversion. 

Figures §(d-f) show the -R v is, cos6* e7 and |cos# 7 | distributions of e7 events obtained using 
the entire data set with cuts applied in the order described above. The irreducible background 
consists almost entirely of e + e~ events with one radiated photon. 

3.4 Kinematic Fits 

The existence of excited leptons would reveal itself as an excess in the total number of observed 
events, appearing as a peak in the reconstructed £7 invariant mass distributions. Kinematic fits 
are used to improve the reconstructed mass resolution of the selected events and also further 
reduce the background thereby increasing the sensitivity of the analysis to excited leptons. 

The kinematic variables used as input to the fit are the energy and direction (E, 9, 0) of 
each identified jet in an event. The energy of tau candidates is left as a free parameter and 
the direction is taken to be the jet axis. A kinematic fit also requires as input the error on 
each measured variable. Estimates of the uncertainties on the energy and direction for the 
different types of leptons and for photons are obtained from studies of di-lepton events in data 
recorded at centre-of-mass energies near and greater than the Z° mass. The error estimates 
are parameterised as functions of the jet energy and polar angle. The uncertainty on the jet 
energy is typically 2 GeV for electrons and photons, and about 5 GeV for muon candidates. The 
uncertainty on the jet polar angle is about 2 mrad for electrons and muons, 4 mrad for photons, 
and 7 mrad for tau candidates. Finally, the azimuthal angle of electron, muon, photon and tau 
candidates is typically known to 0.4 mrad, 0.4 mrad, 3.5 mrad and 7 mrad respectively. 

The kinematic fit enforces conservation of energy and momentum while taking into account 
the beam energy spread as measured by the LEP energy working group ^ . This last constraint 
is necessary since the expected mass resolution for excited leptons is of the same order as the 
centre-of-mass energy spread, which is measured to be approximately 250 MeV. 

Slightly different kinematic fits are applied for single and pair production candidate events. 
In addition to energy and momentum conservation, the kinematic fits for pair produced excited 
lepton candidates also require the invariant masses of the two lepton-photon pairs in the event 
to be equal. There are two possible lepton-photon pairings in each event, and for each pairing 
an additional fit is performed assuming the presence of an undetected initial state radiation 
photon along the beam axis. Thus for each pair production candidate, four kinematic fits are 
performed. Similarly, two kinematic fits are applied to singly produced excited lepton candidate 
events. In the first case, only the two leptons and one photon are included in the fit. In the 
second case, the fit is performed assuming the presence of an initial state radiation photon along 
the beam axis. Finally, a single kinematic fit is performed for e7 events. The fit assumes the 
presence of an undetected electron along the beam axis. 

For a given final state, events are rejected if every kinematic fit attempted has a probability 
less than 0.001. When more than one successful kinematic fit is obtained for an event, the fit per- 
formed without the presence of an initial state radiation photon is chosen if the fit probability is 
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greater than 0.001. For pair production candidates where two fits without initial state radiation 
are performed, the lepton-photon pairing corresponding to the fit with the highest probability 
in excess of 0.001 is chosen. Otherwise, results from the fit with an initial state radiation photon 
are retained where, for pair production candidates, the lepton-photon pairing corresponding to 
the fit with the highest probability is chosen. This requirement on the kinematic fit probabilities 
reduces the number of selected events by more than 70% for ££77 final states and by about 10% 
for final states compatible with the single production of an excited lepton. 

Results from the chosen kinematic fit for each event are used to calculate £7 invariant masses. 
Using the procedure outlined above, the correct lepton-photon pairing for pair produced excited 
leptons is chosen more than 98% of the time as determined using simulated signal events. For £lj 
events, two lepton-photon combinations are possible, both of which are included in the analysis. 
Mass resolutions of approximately 0.2-0.4 GeV for excited electrons and muons and 0.7-2.0 GeV 
for excited taus are obtained using results from the kinematic fits. The natural decay width 
of excited leptons for couplings not excluded by previous searches is constrained to be much 
smaller than these mass resolutions and is thus neglected. Figure || shows the invariant mass 
distributions of selected tt^ and e7 events. There are two entries per li-y event, corresponding 
to the two possible lepton-photon pairings. These distributions are obtained by combining 
data from all the centre-of-mass energies considered. In both single and pair production event 
samples, no mass peaks are observed in the data. 

4 Results 

The numbers of events observed in the data and the corresponding numbers of background events 
expected from Standard Model processes are shown in Table |2|. Selected candidate events for 
the pair production of excited leptons are listed in Table ||[ Typical selection efficiencies for the 
pair production of excited leptons vary from about 35% to 55%. The efficiency for the single 
production of excited muons is 70% and approximately constant over the entire kinematically 
allowed range of masses. Near the kinematic limit for the single production of excited taus, 
the efficiency rapidly drops from 53% down to approximately 20% since the recoiling tau has 
low energy and thus often fails the initial set of selection criteria. For singly produced excited 
electrons, the efficiencies of the ee7 and e7 selections depend strongly on the mixture of s-channel 
and i-channel components. The sum of the ee7 and e7 efficiencies is typically between 50% and 
70%. The evaluation of the systematic uncertainties on the selection efficiencies and background 
estimates is discussed in the following section. No excess of data indicating the existence of 
excited leptons is found in either the single or pair production search. 

4.1 Systematic Uncertainties 

The following sources of systematic uncertainties on the signal efficiencies and background esti- 
mates were investigated. These are described in order of importance. 

Uncertainties in the modelling of initial state photon radiation (ISR) in di-lepton events affect 
the background estimates. They are assessed by comparing background expectations from the 
KORALZ and KK2F event generators for the processes e + e~ — > and e + e~ — > t + t~ . The 

Monte Carlo program KK2F, used in this analysis to estimate the background contributions from 
fi^ 'fi~ and t + t~ events, has the most complete description of initial state radiation including 
second-order subleading corrections and the exact matrix elements for two hard photons [pSfl . 
The relative variations in background expectations between the two Monte Carlo generators are 
assigned as systematic uncertainties representing the effect of missing higher orders. These are 
found to be 11% for final states compatible with the single production of excited muons and 
taus, and 7% for /U//77 and TT77 events. The BHWIDE and TEEGG event generators, used to 
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Table 2: Total numbers of selected events in the data and expected numbers of background 
events for the different final states considered. Statistical and total systematic uncertainties on 
the background estimates are also shown. Contributions to the total systematic errors on the 
background expectations are listed in the lower part of the table. 
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Table 3: List of selected candidate events for the pair production of excited leptons. For each 
candidate, the centre-of-mass energy and reconstructed invariant mass obtained after performing 
the kinematic fits are listed. 



simulate the background from radiative e + e~ events, have a precision for radiative corrections 
similar to the KORALZ program. The background estimates for events expected from the 
production of excited electrons are thus assigned an uncertainty of 7% for the ee77 final state 
and 11% for both ee7 and e7 events. This uncertainty is significantly larger than the error for 
inclusive electron pair production cited in [30|. 

For the purpose of calculating limits on the product of the cross-section and the branching 
fraction, it is necessary to be able to calculate the efficiency and mass resolution of signal events 
at arbitrary excited lepton masses and centre-of-mass energies. For each final state, the selection 
efficiencies and mass resolutions are parameterised as a function of the excited lepton mass 
scaled by the centre-of-mass energy (me* / y/s) . The systematic uncertainties associated with the 
interpolation of efficiencies and mass resolutions were estimated by calculating the root-mean- 
square spread between simulated signal event samples and the parameterisation functions. 

Uncertainties on the fit variable error estimates are evaluated by varying the errors on each 
variable independently. The errors are varied by an amount representing one standard deviation 
as calculated from the uncertainties on the energy and angular parameterisation. Background 
estimates for final states containing two leptons and two photons are particularly sensitive to 
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changes in the errors due to the additional constraint in the kinematic fit requiring events to have 
equal reconstructed lepton-photon invariant masses. Also, the smaller sample of tau pair events 
used to parameterise the errors on the tau direction results in larger statistical uncertainties 
on the error parameterisation which in turn dictate the larger variations used to estimate the 
systematic error contributions. 

The jet classification into leptons or photons contributes to the overall systematic uncertainty 
through the modelling of the lepton and photon identification efficiencies. Using di-lepton and 
di-photon events recorded at centre-of-mass energies equal to and greater than the Z° mass, 
the systematic uncertainty associated with each set of lepton and photon requirements was 
evaluated by comparing the identification efficiencies obtained from data and simulated events. 
Relative errors of 1% for electron and muon, and 2% for the tau and photon classifications are 
assigned. Systematic uncertainties associated with each final state were determined by adding 
linearly contributions from identical jet classifications and adding in quadrature contributions 
from different types of leptons and photons. The resulting uncertainties on the signal efficiencies, 
shown in Table ||, are fully correlated with the corresponding errors on the background estimates 
presented in Table ||. 

The systematic uncertainty associated with the energy scale, energy resolution and angular 
resolution of the leptons and photons was evaluated by modifying each parameter independently 
in Monte Carlo simulated events. Comparisons between data and simulated distributions of 
didepton events recorded at different centre-of-mass energies were used to determine the size 
of these variations. The energy (momentum) of electron and photon (muon) candidates was 
shifted by 0.3%. The energy and angular resolutions of jets were smeared by the maximum 
values for which simulated events were compatible with the distribution of data within one 
standard deviation. Variations in the energy scale result in negligible changes in efficiencies and 
background. Contributions to the systematic uncertainty of each final state from individual 
changes in the energy and angular resolution are added in quadrature. 

The systematic uncertainty due to Monte Carlo modelling of the event selection variables 
was estimated by varying each selection cut independently and measuring the corresponding 
changes in the overall signal efficiencies and background estimates. The difference between the 
mean value of the data and background expectation for each selection variable determined the 
range of variation of each cut. Systematic uncertainties varying between 0.5% and 6.3% are 
assigned to the different background estimates. Contributions to the systematic error on the 
signal efficiencies are shown in Table |||. 

Lastly, the uncertainty on the integrated luminosity measurements (0.2%) is considerably 
smaller than the systematic effects already described and is therefore neglected. 

Summaries of the systematic effects on the background expectations and signal efficiencies 
are presented in Tables [2] and ||, respectively. These systematic uncertainties are included in the 
calculation of limits as described in the following section. 

4.2 Limit Calculations 

Limits on the product of the cross-section and the electromagnetic branching fraction of excited 
leptons are obtained from both pair and single production searches. The numbers of data 
and expected background events at each centre-of-mass energy are binned as a function of the 
reconstructed invariant mass. For selected ££j candidates, both possible i'y invariant masses are 
used. Due to the excellent mass resolution, the double-counting of il'y events does not affect the 
limits calculated. Each mass bin at a given centre-of-mass energy is treated as an independent 
counting experiment. 

For the purpose of calculating limits, the signal invariant mass is assumed to be well described 
by a Gaussian distribution centred at the test mass value and with a width equal to the expected 
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Table 4: Relative systematic uncertainties on the signal efficiencies for each final state consid- 
ered. 

mass resolution. The validity of this assumption is verified with Monte Carlo simulation of 
signal events at different masses and centre-of-mass energies. Efficiency corrections due to non- 
Gaussian tails in the invariant mass distributions are applied to the signal expectation. These 
correction factors, signal efficiencies and mass resolutions are all parameterised as a function of 
the excited lepton mass scaled by the centre-of-mass energy, mi* /y/s. The efficiency correction 
factors are constant over the entire kinematically allowed range and vary from approximately 0.7 
to 0.85 depending on the event final state. Efficiencies and mass resolutions are well-described 
by polynomial functions of various degrees. The efficiencies for the single production of excited 
leptons are calculated with the production and decay angular distributions corresponding to 
/ = /'. The assignment f = f particularly affects the relative fraction of excited electron 
events in the e7 and ee7 selections. 

For a given test mass, the Gaussian distributions describing the invariant mass of signal 
events at each centre-of-mass energy considered are normalised to the expected excited lepton 
cross-section at the highest centre-of-mass energy, thereby taking into account the energy depen- 
dence of the cross-section. A likelihood ratio method [|3l| is used to compute the 95% confidence 
level upper limit on the number of signal events produced in the entire data set (N95). Sys- 
tematic uncertainties on the signal efficiency and background expectation are incorporated by 
fluctuating, over many iterations, the background expectation and signal efficiency according to 
their respective systematic uncertainties. The final limits are determined from the average of all 
the N95 values obtained at each iteration. Systematic errors on the background estimates are 
treated as being fully correlated. The systematic uncertainties on the signal efficiencies due to 
the jet classification are also fully correlated with the corresponding errors on the background 
estimates and are treated as such in the limit calculations. 

Limits on the product of the cross-section and the branching fraction are scaled to y/s = 
208.3 GeV assuming the cross-section evolution as a function of centre-of-mass energy expected 
for excited leptons. The upper limits on the single production of excited muons and tau lep- 
tons do not depend on the coupling assignment of / and /'. The excited electron selection 
efficiencies, however, depend on the relative magnitude of the s-channel and i-channel diagrams. 
For comparison with previously published results, the limits on excited electrons presented here 
assume / = /'. Figures ||(a,b) show the 95% confidence level upper limits on the product of the 
cross-section at y/s = 208.3 GeV and the branching fraction obtained from the search for singly 
and pair produced excited leptons. 

The theoretical calculation H of the product of the pair production cross-section at y/s = 
208.3 GeV and the branching fraction squared is overlayed on Figure |](b). As part of this 
calculation, the electromagnetic branching fraction is calculated assuming / = /'. The 95% 
confidence level lower mass limits on excited leptons correspond to the mass at which the cross- 
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section times branching fraction limit curves cross the theoretical expectation. Lower mass limits 
of m e * > 103.2 GeV, m M » > 103.2 GeV and m T * > 103.2 GeV are obtained. Although systematic 
errors are incorporated into the limit calculations, an additional uncertainty on the mass limits 
arises from the finite width of the centre-of-mass energy bins considered. The 0.5 GeV centre- 
of-mass energy bin width near the kinematic limit corresponds to an uncertainty of 0.1 GeV on 
the mass limits. 

Limits on the product of the cross-section and the electromagnetic branching fraction of 
singly produced excited leptons are used to constrain parameters of the model introduced in 
Section |l[ Since the cross-section for the single production of excited leptons is proportional to 
(//A) 2 , limits on the ratio of the coupling to the compositeness scale as a function of excited 
lepton mass are calculated using 



where N exp is the number of expected signal events assuming //A = 1 TeV" 1 . Figure |(c) shows 
these limits for each type of excited lepton. The //A limit for excited electrons is approximately 
an order of magnitude better than for muons and taus due to the enhancement of the cross- 
section coming from the ^-channel contribution. 

5 Conclusion 

A search for electromagnetically decaying charged excited leptons was performed using 684.4 pb _1 
of data collected by the OPAL detector at \/s = 183-209 GeV. No evidence was found for the 
existence of excited leptons. Upper limits on the product of the cross-section and the branching 
fraction were calculated. From pair production searches, 95% confidence level lower limits on 
the mass of excited leptons are determined to be mi* > 103.2 GeV for t = e, ^,r. From the 
results of the search for singly produced excited leptons, limits were calculated on the ratio of 
the coupling constant to the compositeness scale (//A) as a function of excited lepton mass. 
The results are currently the most stringent constraints on the existence of excited leptons and 
therefore represent a significant improvement on limits previously published |2|-Q. 

Acknowledgements 

We particularly wish to thank the SL Division for the efficient operation of the LEP accel- 
erator at all energies and for their close cooperation with our experimental group. In addition 
to the support staff at our own institutions we are pleased to acknowledge the 
Department of Energy, USA, 
National Science Foundation, USA, 
Particle Physics and Astronomy Research Council, UK, 
Natural Sciences and Engineering Research Council, Canada, 

Israel Science Foundation, administered by the Israel Academy of Science and Humanities, 
Benoziyo Center for High Energy Physics, 

Japanese Ministry of Education, Culture, Sports, Science and Technology (MEXT) and a grant 

under the MEXT International Science Research Program, 

Japanese Society for the Promotion of Science (JSPS), 

German Israeli Bi-national Science Foundation (GIF), 

Bundesministerium fur Bildung und Forschung, Germany, 

National Research Council of Canada, 




12 



Hungarian Foundation for Scientific Research, OTKA T-029328, and T-038240, 
Fund for Scientific Research, Flanders, F.W.O.-Vlaanderen, Belgium. 



References 

[f] F.E. Low, Phys. Rev. Lett. 14 (1965) 238. 

[2] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. C14 (2000) 73; 
OPAL Collaboration, K. Ackerstaff et al, Eur. Phys. J. CI (1998) 45; 
OPAL Collaboration, K. Ackerstaff et al, Phys. Lett. B391 (1997) 197; 
OPAL Collaboration, G. Alexander et al, Phys. Lett. B386 (1996) 463; 
OPAL Collaboration, M.Z. Akrawy et al, Phys. Lett. B244 (1990) 135. 

[3] L3 Collaboration, M. Acciarri et al, Phys. Lett. B502 (2001) 37; 
DELPHI Collaboration, P. Abreu et al, Eur. Phys. J. C8 (1999) 41; 
ALEPH Collaboration, R. Barate et al, Eur. Phys. J. C4 (1998) 571. 

[4] HI Collaboration, C. Adloff et al, Phys. Lett. B525 (2002) 9; 
HI Collaboration, C. Adloff et al, Eur. Phys. J. C17 (2000) 567; 



ZEUS Collaboration, S. Chekanov, et al, Submitted to Phys. Lett. B [ hep-ex/0109018 ] 
ZEUS Collaboration, J. Breitweg et al, Z. Phys. C76 (1997) 631; 
ZEUS Collaboration, M. Derrick et al, Z. Phys. C65 (1995) 627. 

[5] F. Boudjema, A. Djouadi and J.L. Kneur, Z. Phys. C57 (1993) 425. 

[6] K. Hagiwara, D. Zeppenfeld and S. Komamiya, Z. Phys. C29 (1985) 115. 

[7] OPAL Collaboration, S. Anderson et al, Nucl. Instr. Meth. A403 (1998) 326; 
B.E. Anderson et al, IEEE Trans. Nucl. Sci. 41 (1994) 845; 
OPAL Collaboration, K. Ahmet et al, Nucl. Instr. Meth. A305 (1991) 275. 

[8] LEP Energy Working Group, A. Blondel et al, Eur. Phys. J. Cll (1999) 573. 

[9] R. Tafirout and G. Azuelos, Comp. Phys. Comm. 126 (2000) 244. 

[10] A. Djouadi, Z. Phys. C63 (1994) 317. 

[11] S. Jadach, W. Placzek and B.F.L. Ward, Phys. Lett. B390 (1997) 298. 

[12] D. Karlen, Nucl. Phys. B289 (1987) 23. 

[13] S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 79 (1994) 503. 

[14] S. Jadach, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 130 (2000) 260; 
S. Jadach, B.F. Ward and Z. Was, Phys. Lett. B449 (1999) 97. 

[15] T. Sjostrand, Comp. Phys. Comm. 82 (1994) 74. 

[16] S. Jadach, W. Placzek, M. Skrzypek, B.F.L. Ward and Z. Was, Comp. Phys. Comm. 
(1999) 272. 

[17] J. Fujimoto et al, Comp. Phys. Comm. 100 (1996) 128. 
[18] F.A. Berends, R. Kleiss, Nucl. Phys. B186 (1981) 22. 



13 



[19] R. Bhattacharya, J. Smith and G. Grammer, Phys. Rev. D15 (1977) 3267; 
J. Smith, J.A.M. Vermaseren and G. Grammer, Phys. Rev. D15 (1977) 3280. 

[20] R. Engel and J. Ranft, Phys. Rev. D54 (1996) 4244. 

[21] G. Marchesini et al, Comp. Phys. Comm. 67 (1992) 465. 

[22] J. Allison et al, Nucl. Instr. Meth. A317 (1992) 47. 

[23] OPAL Collaboration, G. Alexander et al.., Z. Phys. C72 (1996) 191. 

[24] OPAL Collaboration, G. Alexander et al, Z. Phys. C52 (1991) 175. 

[25] OPAL Collaboration, R. Akers et al, Z. Phys. C63 (1994) 197. 

[26] OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. C16 (2000) 41; 
OPAL Collaboration, G. Abbiendi et al, Eur. Phys. J. C8 (1999) 217; 
OPAL Collaboration, G. Alexander et al, Z. Phys. C70 (1996) 357. 

[27] OPAL Collaboration, K. Ackerstaff et al, Phys. Lett. B389 (1996) 416. 

[28] OPAL Collaboration, K. Ackerstaff et al, Eur. Phys. J. C2 (1998) 213. 

[29] S. Jadach, B. F. L. Ward and Z. Was, Phys. Rev. D63 (2001) 113009. 

[30] S. Jadach, W. Placzek and B.F.L. Ward in "Two-Fermion Production in Electron-Positron 
Collisions", edited by M. Kobel and Z. Was, CERN-2000-009-D jhep-ph/000718C|l . 

[31] T. Junk, Nucl. Instr. Meth. A434 (1999) 435. 



14 



OPAL 



3 i 

| 10 3 

<*> in 2 

.2 1" 

a 10 
- 

1 



10 
10 



(a) eeyy 



i i i i I i i i i I i i i i I i i i i i i i i 
(d) eey ^ 




0.25 0.5 0.75 1 



R 



0.25 0.5 0.75 1 



R 



3 i 
| 10 3 

I io 2 

1 10 

- 

1 



10 
10 



(b) w^yy r> 



' ' ' ITVHi I 'i I r 



0.25 0.5 0.75 




- 1 1 1 1 1 1 1 1 1 1 1 

r (e) |W 




1 1 1 1 1 1 1 - 

r> i 


-+■*■ 

_ f 

1 T 

■< 




-If!*!*' 

f t 


r ' . : ; ; r 


I i i i 


1 1 1 1 1 1 1 




: i 1 1 1 i i'v 



0.25 0.5 0.75 1 



R 



R 



-r 

© i 
^ 1(1 

.2 10 1 

I 10 

1 



(c)xxyy 



10 
10 



-1 




I 10 3 

& 

W 10 



0.25 0.5 0.75 1 



10 



-l 



- 1 1 1 1 1 1 1 

r (f) xxy 
pi* 




1 | 1 1 1 - 






*% 


i i i i i 'i i 


1 1 1 1 1 1 1 1 1 


i 1 : i u_ 1 



0.25 0.5 0.75 1 



R 



R 



Figure 1: Distributions of the sum of the energies of the two leptons and one or two photons 
divided by the centre-of-mass energy for pair (a-c) and single (d-f ) production candidate events 
after the preselection. The points represent the combined data from all centre-of-mass ener- 
gies considered while the solid lines are the total expected background from Standard Model 
processes. The dashed lines represent an example of 1*1* (a-c) and £*(d-i) signal events with 
arbitrarily chosen masses of 40 GeV and 90 GeV respectively. The dotted lines show the ex- 
pected distributions for pair and singly produced excited leptons with masses of 90 GeV and 
180 GeV, respectively. Distributions of excited lepton signal events are normalised to a ratio of 
the coupling constant to the compositeness scale of 1 TeV -1 . The arrows indicate the accepted 



regions. The background modelling for low values of R v i s is discussed in Section I3T 
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Figure 2: Distributions of selection variables for different final states: (a) cosine of the angle 
between the most energetic electron and the photon in ee7 events, (b) total energy deposited in 
the electromagnetic calorimeter scaled by the centre-of-mass energy for TT7 events, (c) absolute 
value of the cosine of the missing momentum vector polar angle for TT7 events, (d) sum of the 
energy of the electron and photon divided by the centre-of-mass energy for e7 events, (e) cosine 
of the angle between the electron and photon for e7 events and (f ) absolute value of the cosine 
of the photon polar angle in e7 events. All selection cuts have been applied, in the same order as 
described in the text, up to that on the variable plotted. The dashed and dotted lines represent 
examples of excited lepton signal events with arbitrarily chosen masses of 90 GeV and 180 GeV, 
respectively, and branching fraction calculated assuming a ratio of the coupling constant to the 
compositeness scale of 1 TeV -1 . 



16 



OPAL 




> 

o 

o 



a 



100 



80 



60 



40 



20 




50 100 150 200 



50 100 150 200 



m eY (GeV) 



m eY (GeV) 




m^(GeV) 



m Ty (GeV) 



Figure 3: Reconstructed invariant mass distributions for (a) ee7, (b) &y, (c) /z/ry and (d) ttj 
candidates after all cuts are applied. The points are data and the solid lines represent the total 
expected background from Standard Model processes. The shaded histograms represent excited 
lepton signal events with an arbitrarily chosen mass of 150 GeV and normalised to a ratio of 
the coupling constant to the compositeness scale of 0.4 TeV -1 (a,b) and 2 TeV -1 (c,d). There 
are two entries per event in (a,c,d) corresponding to the two possible pairings. 
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Figure 4: The 95% confidence level upper limits on the product of the cross-section at 
y/s = 208.3 GeV and the branching fraction for (a) single and (b) pair production of excited 
leptons as a function of mass (m*). The limit obtained for the single production of excited elec- 
trons is calculated assuming / = /'. The regions above the curves are excluded. The product of 
the theoretical cross-section and the branching fraction squared assuming / = /' is also shown in 
(b). The 95% confidence level upper limits on the ratio of the excited lepton coupling constant 
to the compositeness scale, //A, as a function of the excited lepton mass and assuming f = f 
are shown in (c). The regions above the curves are excluded by single production searches while 
pair production searches exclude masses below 103.2 GeV for excited electrons, muons and taus. 



18 



OPAL 




100 120 140 160 180 200 220 



m ;: (GeV) 



OPAL 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 


i i I i i 


3 10 










- 








CQ 






X 






1 


e e 

* * 






* * 






X X 




-l 

10 


Theory (f=f) 1 




-2 


L. ZJ 




10 


i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 


i i i i i 



40 50 60 70 80 90 100 110 



m ;: (GeV) 



OPAL 




100 120 140 160 180 200 220 



m ;: (GeV) 



